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Abstract. Supersymmetry with breaking of R-parity provides an attractive way to 
generate neutrino masses and lepton mixing angles in accordance to present neutrino 
data. We review the main theoretical features of the bilinear R-parity breaking (BRpV) 
model, and stress that it is the simplest extension of the minimal supersymmetric 
standard model (MSSM) which includes lepton number violation. We describe how 
it leads to a successful phenomenological model with hierarchical neutrino masses. 
In contrast to seesaw models, the BRpV model can be probed at future collider 
experiments, like the Large Hadron Collider or the Next Linear Collider, since the 
decay pattern of the lightest supersymmetric particle provides a direct connection 
with the lepton mixing angles determined by neutrino experiments. 
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1. Introduction 



A combination of solar, atmospheric, reactor and accelerator neutrino experiments 0, 
0, 0] have now firmly established the existence of neutrino masses and therefore the 
incompleteness of the standard model of electroweak interactions. The determination of 
neutrino oscillation parameters presented in Ref. uses the most recent data and state- 
of-the-art solar and atmospheric neutrino fluxes. For previous reviews and references 
jE S We have now learned that the atmospheric oscillations involving <-» v T 



sec 

are characterized by a nearly maximal mixing, while the solar neutrino mixing angle is 
large, but significantly non-maximal. With the recent standard solar model fluxes there 
is a unique range for the solar mass splitting Am^, determined from the data to be 
about 30 times smaller than the atmospheric mass splitting Am| TM . 

The discovery of neutrino mass constitutes the only solid hint we currently have 
of physics beyond the standard model. There are theoretical arguments based on the 
stability of the gauge hierarchy which suggest the existence of physics at the TeV scale. 
Supersymmetry una provides an answer to both these issues which fits well with 



unification and string theory ideas 11 



Prompted by these data there has been a rush of theoretical and phenomenological 
papers on models of neutrino masses and mixings. The most popular idea is to ascribe 
neutrino masses to physics ^at a large mass scale in order to implement some variant 

13, Ej, 13 • Broken R-parity supersymmetry provides 



of the see-saw mechanism [12, 



a theoretically interesting and phe nomenologically viable alternative to the origin of 



neutrino mass and mixing la, LL2 



221. Here we focus on the case of 



supersymmetry with bilinear R-parity breaking |23j . This is the simplest of all R parity 
violating models. It also provides the simplest extension of the MSSM to include 
the violation of lepton number, as well as a calculable framework for neutrino masses 
and mixing angles in agreement with the experimental data In this 

model the atmospheric neutrino mass scale is generated at the tree-level, through the 
mixing of the three neutrinos with the neutralinos, in an effective 'low-scale" variant 
of the seesaw mechanism 0. la contrast, the solar mass and mixings are generated 
radiatively 26] . BRpV can be considered either as a minimal extension of the MSSM 
|2sL Eoi 0, (with no new particles) valid up to some very high unification energy 
scale, or as the effective description of a more fundamental theory in which the breaking 
of R-parity occurs in a spontaneous way by minimizing the scalar potential 

This short review is mainly devoted to the generation of neutrino masses and lepton 
mixing, both the tree-level atmospheric neutrino mass scale as well as a description of 
the main features of the full one-loop calculation of the neutrino-neutralino mass matrix 
and its various analytic approximations which, in some cases, can be rather simple. For 
definiteness we will stick to the case of explicit BRpV only. 

However, in contrast to the seesaw mechanism, in broken R-parity supersymmetry 
neutrino masses are generated at the electro- weak scale 2J, |26J, |2j|. Such low-scale 



schemes for neutrino masses have the advantage of being testable also outside the 
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realm of neutrino experiments. Although neutrino properties can not be predicted 
from first principles, their fit to the data allows for unambi guo us tests of the theory at 

accelerator experiments |27 . EH H H H H H H13- 

Indeed, the measured 



lepton mixing angles lead to well defined predictions of the decay properties of the 
lightest supersymmetric particle (LSP). This is a very general and robust feature of 
these theories, which holds irrespective of the nature of the LSP. Here we will illustrate 
possible phenomenological scenarios by discussing some examples of measurements of 
decay properties of different LSP candidates. 

This paper is organized as follows. In Sec. El we introduce the main features of the 
model, discuss the soft supersymmetry breaking terms, as well as the relevant fermion 
mass matrices and the main features of the corresponding diagonalizing matrices. In 
Sec. EH we discuss the generation of the atmospheric neutrino mass scale at the tree-level, 
while in Sec. HI we analyse the main features of the one-loop-induced solar neutrino mass 
scale, including a discussion of the relevant Feynman graph topologies. We also give 
simplified approximation formula for the solar mixing angle. We then turn briefly to 
collider phenomenology and how the model under discussion could be tested in LSP 
decays in Sec. El before we conclude and summarize our results in Sec. El 

2. Formalism 

In this section we introduce the main features of the model and the relevant mass 
matrices. The superpotentiai of the model and the soft SUSY breaking terms are given, 
approximate solutions to the tadpole equations discussed. 

2.1. The Superpotentiai and the Soft Breaking Terms 

The minimal BRpV model we are working with is characterized by the presence of three 
extra bilinear terms in the superpotentiai analogous to the \x term present in the MSSM. 
Using the conventions of ref. [3jJ it may be given as 

W = e ab [h%Q a i U j H h u + h%Q\D 3 H a d + h%L%H a d - f,H a d H b u + e&H*] (1) 

where the first three terms are the usual MSSM Yukawa terms, [i is the Higgsino mass 
term of the MSSM, and e« are the three new terms which violate lepton number in 
addition to R-Parity. The couplings hu, ho and hs are 3x3 Yukawa matrices and \i and 
Cj are parameters with units of mass. The smallness of the bilinear term e , in eq. (0) may 
arise from a suitable symmetry. In fact any solution to the jx problem [44] potentially 
explains also the "ej-problem" 0|. A common origin for the 6j terms that account 
for the neutrino oscillation data, and the \x term responsible for electroweak symmetry 
brea king can be ascribed to a horizontal family symmetry of the type suggested in 
Ref. jjj- 

The smallness of e; could also arise dynamically in models with spontaneous 



breaking of R parity |32j, |33|, |34j , where it is given as the product of a Yukawa coupling 



times a singlet sneutrino vacuum expectation value. 
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Supersymmetry breaking is parameterized with a set of soft supersymmetry 
breaking terms. In the MSSM these are given by 

c mssm = MfQTQ* + MfUJJ* + MfDiD* + Mfl^L] + Nlf^R* 
+ m 2 H H a d *H a d + m 2 Hu H a u *H a u - [\M S X S \ S + ±MAA + \M'\'X + h 



d 

^Q-U 3 H b u + A%Q\D 3 H a d + AtLfaHS - BpHSHl 

MSi 
soft 



In addition to the MSSM soft SUSY breaking terms in £^ff M the BRpV model contains 
the following extra terms 

v s B oft V = -Bi€ie ab LfHl , (3) 

where the Bi have units of mass. In what follows, we neglect intergenerational mixing 
in the soft terms in eq. (J2J). 

The electroweak symmetry is broken when the two Higgs doublets H d and H u , and 
the neutral component of the slepton doublets L\ acquire non-zero vacuum expectation 
values (vevs). These are calculated via the minimization of the effective potential or, in 
the diagrammatic method, via the tadpole equations. The full scalar potential at tree 
level is 



v° = V 

' total / j 



dW 2 



\T , V MSSM , TrBRpV r a\ 

V D + V S oft + V soft I 4 J 



where is any one of the scalar fields in the superpotential in eq. (JTJ), Vd are the 
D-terms, and V go ff is given in eq. ©• 

The tree level scalar potential contains the following linear terms 

VLear = M + « + + « + , (5) 

where the different t° are the tadpoles at tree level. They are given by 
t°d = ( m L + l*\vd + v d D - /i(Bv u + 
t° u = - Bfi,v d + (m 2 Hu + fi 2 ^jv u - v u D + ViB^ + v u e 2 
t\ = v x D + e x ( - fiv d + v u Bx + v l e?j + \ (v t M 2 Lil + M 2 u v^\ (6) 
t% = v 2 D + e 2 (- fiv d + v u B 2 + u^) + \ [v t M\ l2 + M 2 L2iVi 
t° 3 = v 3 D + e 3 ( - /iv d + v u B 3 + u^) + \ (y t M 2 l3 + M 2 3iVi 

where we have introduced the notation 

-®- 

and shifted the neutral fields with non-zero vevs as 

H° d = -L[a° d + Vd + i(Pd } , Hi = -L[a u +v u +i<p° u } , Z° = -L[^+v i+ i^} .(8) 
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The five vacuum expectation values can be expressed in spherical coordinates as 

v d = v sin 6*i sin 9 2 sin 9 3 cos (3 
v u = v sin 8i sin 9 2 sin 9 3 sin j3 

v 3 = v sin 9\ sin 9 2 cos 9 3 (9) 
v 2 = f sin 61 cos #2 
ui = v cos 6*i 

which preserves the MSSM definition tan/5 = v u jv& with the boson mass given as 
\g 2 {vj+vl+v 2 +v 2 +v 3 ). We have also defined D = Ug 2 +g' 2 ){vf+v 2 +vl+vj—vl) 



m 



w — 4 

and e 2 = ef+e 2 +e 3 . A repeated index i in eq. (J0J) implies summation over i = 1,2, 3. The 
five tree level tadpoles t° are equal to zero at the minimum of the tree level potential, 
and from there one can determine the tree level vacuum expectation values. 

2.2. Radiative Breaking of the Electroweak Symmetry 

A reliable description of electroweak symmetry breaking and Higgs boson physics in 
supersymmetry requires the inclusion of radiative corrections. In the BRpV model the 
full scalar potential at one-loop level, called effective potential, is 

Vtotal = V t ° otal + V RC (10) 

where K rt n f is given in Eq. (0J and Vrc include the quantum corrections. Following 
Refs. ' 



241 . |26| we use the diagrammatic method, incorporating the radiative corrections 
through the one-loop corrected tadpole equations. The one loop tadpoles are 

t a = tl- St™ + T a (Q) = t + T™{Q) (11) 

where a = d,u, 1,2,3 and T® R (Q) = —5t^ IS + T a (Q) are the finite one loop tadpoles. 
At the minimum of the potential we have t a — 0, and the vevs calculated from these 
equations are the renormalized vevs. 

Neglecting intergenerational mixing in the soft masses, the five tadpole equations 
can be conveniently written in matrix form as 

K 4 *?> 4 4] T = M t 2 ad K v d , Vl ,v 2 , v 3 f (12) 



where the matrix M^ ad is given in 26( and depends on the vevs only through the D 
term defined above. 

In the MSSM limit, where = Vi = 0, the angles 9i are equal to 7r/2. In addition 
to the above MSSM parameters, our model contains nine new parameters, e^, Vi and 
Bi. Considering we have three tadpole equations one can take either the 3 B^ as input 
and derive the 3 sneutrino vevs or vice versa, such that we have in total just six new 
parameters (compared to the MSSM). 

In order to have approximate solutions for the tree level vevs, consider the following 
rotation among the and lepton superfields: 

M' t 2 ad = RM 2 ad R 1 (13) 
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R 



where the three angles are defined as 
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c l = — 



C'2 



C3 



6l 



£2_ 

e 3 



S2 

s 3 



fi 2 + e 



i ■ 



/^ ,2 + 4 , 



/r" 



/i" 2 + el 



(15) 



It is clear that this rotation R leaves the D term invariant. The rotated vevs are given 
by 



K> <4 V 'n V 2i V 3] T = R K> ^ u 2, W 3] T 



(16) 



and under the assumption that v[,v' 2 ,v' 3 <C v, these three small vevs have the 
approximate solution 

mi. - M? 



M' 2 +D 



n'n'" 



Bx-B , 
-v' 



//e 2 



M' 2 +D 



m 



fJL"fJL'" 



Mt 2 , B 2 



B' , 

— v„ 



Mg + D 



in 



1/2 



-M\ 



^112 p„ 

where we have defined the following rotated soft terms: 



B, - B" 



(17) 



™ 2 H y + M 2 Li e 2 



B' 



M 



r2 



/i' 2 

By 2 + B x e\ 

n' 2 

™\A + M 2 y 



m 



112 



m' 2 df i' 2 + M 2 2 e 2 



112 



m 



12 



B" 



B'fx' 2 + B 2 e 2 



ml 



m" 2 y 2 + M 2 3 e 2 



1/12 



112 



m' 2 d + M[ a 



12 



/// 



112 



B" 



M' 



^ + M, (18) 



A 4 ' 



1112 



nr 



'SA + M IA' 



1112 



As can be seen from eq. (17) the approximation v'^v^v'^ <C v is justified if either 
a) £i < and/or b) (m 2 Hd - M 2 Li )j^ 2 < 1 and (5, - B)/fi < 1. The latter holds 
automatically (to some extent) in many models of supersymmetry breaking, as for 
example in minimal supergravity Ref. [23]. 

As in the MSSM, the electroweak symmetry is broken because the large value of the 
top quark mass drives the Higgs mass parameter m 2 Hu to negative values at the weak 
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scale via its RGE 47]. In the rotated basis, the parameter fi'" 2 is determined at one 
loop by 



1112 



rrir 



A zz (m%) 



+ 



JII2 



DR 



m 



DR 



t 



1 



(19) 



where t'a = v' u /v' d is defined in the rotated basis and is analogous to tan/3 in eq. © 
defined in the original basis. The finite DR Z-boson self energy is A^^(m|), and the 
one-loop tadpoles T® R and T® R are obtained by applying to the original tadpoles in 
eq. the rotation R defined in eq. (jl4J) . The radiative breaking of the electroweak 
symmetry is valid in the BRpV model in the usual way: the large value of the top quark 
Yukawa coupling drives the parameter m 2 Hu to negative values, breaking the symmetry 
of the scalar potential. 



2.3. Neutral fermion mass matrix 

Here we consider the tree level structure of the fermion mass matrices in this model. 
For a more complete discussion of different mass matrices in BRpV see the Appendix 
of Ref. 0- In the basis ^ 0T = HA' 
matrix Mjv is given by 



i\ 3 , Hj, H 2 , z/ e , z/ T ) the neutral fermion mass 



M 



N 



m 



(20) 



where 



Mi 




- \g'v d 



M 2 

\gv d 



-\g'v d \g'v u 
\gvd -\gv u 

-fi 



(21) 



-/i 

is the standard MSSM neutralino mass matrix (M 2 and Mi are the SU(2) and U(l) 
gaugino soft masses) and 

l 



m 



ig'vi 



(22) 



iffvi ei 
\g'v 2 \gv 2 e 2 
\g'vz \gv?, e 3 

characterizes the breaking of R-parity. The full 7x7 neutrino/neutralino mass matrix 
is diagonalized as 

N*Wl N M~ l = diag(m x o, m v .) (23) 
where (i — 1, • • • , 4) for the neutralinos, and (j = 1, • ■ ■ , 3) for the neutrinos. 

A^M^A/"- 1 = M d Jo 9 (24) 
and the eigenvectors are given by 

F» = Ni^ (25) 
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using the basis ip 
d ap 

TV* 



—i\ 3 , Hj, H%, u e , v T ). As discussed in more detail below, 



to a very good approximation, the rotation matrix can be written as 

N* JV-ft 



-V T £ V 

y v S " i 



T 



(26) 



Here, N is the rotation matrix that diagonalizes the 4x4 MSSM neutralino mass matrix, 
V v is the rotation matrix that diagonalizes the tree level neutrino 3x3 mass matrix, and 
£ij <C 1 are the relevant small expansion parameters which characterize the violation of 
R-parity and whose form will be given in Sec. 01 

2.4- Charged fermion mass matrix 

The chargino/lepton mass matrix is given by 



M 



c 



M 2 


T2 gVu 













/' 




-71(^)22 ^2 


-75(^)33 ^3 


. 73^ 3 


-ei 
-ea 
-e 3 








72 (M22 w d 







71 (M33 v d 



(27) 



We note that the chargino sector decouples from the lepton sector in the limit e« = Vi = 0. 
As in the MSSM, the chargino mass matrix is diagonalized by two rotation matrices U 
and V defined by 

rdiag 



W M^V" 1 = M 



F+ 



with the eigenvectors satisfying 

-iX + , R\, e%, fi%, r£) and ip- 



F + 

r Ri 



(28) 



(29) 



in the basis ip + 
Dirac fermions being 



-iX , Hf, e L , fi L , r L ), and with the 



/ Ft \ 



\ F ^ J 



(30) 



To first order in the R-Parity violating parameters we have 



V 



v ve R 



u 



R 



u 

-VttL 



vt 



(31) 



where V[* and Vr diagonalize the charged lepton mass matrix according to V[*1S/L £ V R = 
Mj ia9 . For most purposes it is sufficient to take £r = 2X 3, since it is smaller than £ L 
typically by a factor of mi/m SUSY - Note, that we can choose V[* = Vj 



R 



L 3x3 



then have 



til 



a\ Aj + bti 



. We 



(32) 
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and 

l 9 l Q 2v u , 00 \ 

01 = 71X7' a > = -^K~ + (33) 

where A + is the determinant of the 2x2 chargino mass matrix and 

Ai = \iVi + v d ei oc v • (34) 
are the alignment parameters. 

3. Tree— level neutrino mass: the atmospheric scale 

The tree-level contribution to neutrino masses from broken R parity supersymmetry has 
a long history ji^ . Thanks to the Super-K findings 0] we will be interested only in the 
case where the neutrino mass which is determined at the tree level is small, in order 
to account for the atmospheric neutrino data. The above form for is especially 
convenient in this case in order to provide an approximate analytical discussion valid in 
the limit of small 1/ip violation parameters. Indeed in this case we perform a perturbative 
diagonalization of the neutral mass matrix, defining 

£ = m • M-J (35) 

Since the effective RpV parameters are smaller than the weak scale, we can work in a 
perturbative expansion defined by £ <C 1, where £ denotes a 3 x 4 matrix given as 

g'M 2f i 
KlX 2det(M x o) 1 

gM lf i 
^ 2 2det{M x o) 1 

€j {g 2 M 1 +g' 2 M 2 )v 
^ 3 /i + Adet(M x o) 1 
(g*Mi + g' 2 M 2 )v d 

U = ~ 4det(M x o) Al (36) 
From Eq. (j3EJ) and Eq. one can see that £ = in the MSSM limit where 6j = 0, 
Vi = 0. 

If the elements of this matrix satisfy 

« 1 (37) 
then one can use it as expansion parameter in order to find an approximate solution for 
the mixing matrix TV. 

In leading order in £ the mixing matrix M is given by, 

The second matrix above block-diagonalizes the mass matrix M^r approximately to the 
form diag(A^ x o, m e //), where 

m eff = - m ■ M~lm T 
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( A 2 



M l9 2 + M 2 g 
Adet(M x o) 



AJV, 



A e A„ 
A 2 



A e A r \ 
A M A r 



\A e A T A M A 7 

The sub-matrices N and V v diagonalize -M x o and m e ff 
N*M x oN f = diag(m x o), 

V, T m eff V u = diag(0, 0,m u ), 



where 



m,. 



Trim, 



M ig 2 + M 2 g 



/2 



eff) 



■|Ap 



(39) 

(40) 
(41) 

(42) 



Adet(M x o) 

The special form of the neutralino/neutrino mass matrix implies that the effective 
neutrino mass matrix m e f / generated after diagonalizing out the heavy neutralinos has a 
projective form, a feature common to many spontaneous R-parity violating models 0]. 
This implies that only one neutrino acquires a tree level mass, the other two remaining 



massless |4j|. As a result at the tree approximation one can rotate away one of the 
three angles in the matrix V v , leading to 



( 1 

cos #23 - sin # 2 3 
y sin #23 cos6> 23 



X 



^ cos #13 — sin 613 

1 
i sin# 13 cos #13 



(43) 



where the mixing angles can be expressed in terms of the alignment vector A as follows: 

A R 



tan 9 



13 



tan # 23 = — r^. 



(A 2 +A 2 ), 

A, 
A T 

The non-zero tree-level eigenvalue of the neutrino mass matrix is identified with the 
atmospheric mass scale. The calculated AmJ (m can be expressed as a function of the 



(44) 
(45) 




Figure 1. Am^ tm versus the BRpV alignment parameters 



alignment parameter A (left in Fig. [!}, or as function of | A|/ ( -y/ M^y) (right in Fig. [TJ, 
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all of these expressed in GeV. The figure shows that Eq. can be used to fix the 
relative size of R-parity breaking parameters to obtain the correct Am 2 atm . On the other 
hand, as shown in Fig Elthe atmospheric angle can be expressed in terms of A^/A r . Its 
maximality is obtained for ~ A T if A e is smaller than the other two. Let us stress 





|A M /A T | 



Figure 2. The atmospheric angle versus the ratio of BRpV parameters |A M /A T | 



once again that there is no solar mass splitting in the tree approximation so that, as a 
result, the "solar angle" is not defined, as it can be rotated away by redefining the two 
degenerate neutrinos 



141. 



4. One— loop— induced neutrino mass: the solar scale 



As we just saw in the BRpV model the atmospheric mass scale and mixing arises at 
the tree-level. We now discuss the determination of solar neutrino masses and mixings, 
which are both generated radiatively. One-loo p co rrections to the neutrino masses can 
be calculated numerically [26] or analytically [24i ]. While the numerical approach can 
give "exact" results (exact in the sense of being correct up to higher order effects), the 
analytic approach, while being less accurate, gives a better understanding about which 
parameters control the loops and thus the solar neutrino mass and angles in our model. 
The discussion will therefore mainly concentrate on the analytical calculations. 

In principle, in order to find the correct neutrino mixing angles one has to 
diagonalize the one-loop corrected neutralino/neutrino mass matrix. We define 

MP° le = Mjf (Q) + AMjj (46) 

where Mjj R (Q) is the tree-level pole mass and DR indicates the dimensional reduction 
scheme we used in the numerical calculation. One-loop corrections are 



AM 



ngK 2 ) 



+ ng(mj) 



m v oE^(m-) + m v oS 



(47) 



where the symmetrization is necessary to achieve gauge invariance and consistency with 
the Pauli principle. Here U v and H v are the renormalized self-energies. They contain 
products of couplings and the usual Passarino-Veltman functions [50j. 
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Diagonalizing the tree-level neutrino mass matrix first and adding then the one-loop 
corrections before re-diagonalization one finds that the resulting neutrino/neutralino 
mass matrix has non-zero entries in the neutrino/neutrino, the neutrino/neutralino and 
in the neutralino/neutralino sectors. We have found 24] that the most important part 
of the one- loop neutrino masses derives from the neutrino/neutrino sector and that the 
one- loop induced neutrino/neutralino mixing is usually negligible. 

The relevant topologies for the one loop calculation of neutrino masses are then 
illustrated in FigEl Here our conventions are as follows: open circles with a cross 






Figure 3. Topologies for neutrino self-energies in the BRpV supersymmetric model 



inside indicate genuine mass insertions which flip chirality. On the other hand open 
circles without a cross correspond to small R-Parity violating projections, indicating 
how much of an Rp-even/odd mass eigenstate is present in a given Rp-odd/even weak 
eigenstate. In the actual numerical calculation these projections really belong to the 
coupling matrices attached to the vertices. However, given the smallness of Rp- violating 
effects, the pre-diagonalization "insertion-method" proves to be a rather useful tool to 
develop an analytical perturbative expansion and to acquire a simple understanding of 
the results. 

These topologies have then to be "filled" with all relevant combinations of 
particles/sparticles. Here we will concentrate on discussing the loop involving 
bottom/sbottom quarks, since a) it is in large parts of parameter space the numeri cally 
most important one. And b) other loops can be calculated in a very similar way [24], 
although they are more complicated. 

The relevant Feynman rules for the bottom-sbottom loops are, in the case of left 
sbottoms: 

bj / _ ■ r/nbn6 (l-7s) i fjmb (l+7s) l 

— l Y u Lij 2 ^ U Rij 2 J 
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with 

(M = - - lLp= tan 9 W ^ 



O h i 3 = R% -jL (M l2 - | tan e w Nn) - i? 2 M? 3 (48) 

After approximating the rotation matrix M we find that expressions with the 
replacement N — > N are valid when the neutral fermion is a neutralino. When the 
neutral fermion F° is a neutrino, the following expressions hold 

of b « (a 3 |A|^ 3 + fee*/) + i^ 2 ^tan^oi|X|^ 3 

°Rvi « ^ii^f (itan^O! - a 2 ) |A|^ 3 + i^/fc (a 3 |A*|fc 3 + (49) 

where z' = i— 4 label one of the neutrinos. _Rj fc are the rotation matrices connecting weak 
and mass eigenstate basis for the scalar bottom quarks. In case of no intergenerational 
mixing in the squark sector R b - k can be parameterized by just one diagonalizing angle 

(>,,■ 

Putting these couplings together one finds the simplest contribution to the 
radiatively induced neutrino mass from loops involving bottom quarks and squarks [3] 

5«(0) = E + 0^ r O b £) m b B (0, ml m 2 r ) (50) 



where B (0,ml,m 2 ) is a Passarino-Veltman function [50j can be written as follows 

N c rr, 

167T 



% = " ^ 2s nh 2 ABl^ x (51) 



+ a 3 b (e^ 3 + i 3 5 i3 ) |A| + (a* + S i3 6 j3 \A\' 

This expression is proportional to the difference of two B functions, 

AB~f b * = B (0,m 2 b ,ml) - B o {0, m 2 , mfj (52) 

Parameters Aj have been defined above. The e parameters are defined as q = (V^Y 3 ej, 
and are given by 

e e (A 2 + A 2 T ) - A e (A^ + A T e T ) 



A r e M - A M e T 



= (53) 



A ■ e 

^3 



+ AJ + A? 
On the other hand a^R are defined as 

a R = Q%cti - a 2 J , a L = |%ai (54) 

The different terms in eq. f|51|) can be understood as coming from the graphs 
corresponding to the first topology of Fig. El They have been depicted in more detail in 
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Figure 4. Bottom-Sbottom diagrams for solar neutrino mass in the BRpV model 



Fig. HI where we have adopted the following conventions: a) as before, open circles 
correspond to small R-parity violating projections, indicating how much of a weak 
eigenstate is present in a given mass eigenstate, (b) full circles correspond to R-parity 
conserving projections and (c) open circles with a cross inside indicate genuine mass 
insertions which flip chirality. 

The open and full circles should really appear at the vertices since the particles 
propagating in the loop are the mass eigenstates. We have however separated them to 
better identify the origin of the various terms. There is another set of graphs analogous 
to the previous ones which corresponds to the heavy sbottom. They are obtained from 
the previous graphs making the replacement bi — > b 2 , s^ — > c~ h and c~ h — > — sy. Note 
that for all contributions to the 2x2 sub-matrix corresponding to the light neutrinos 
the divergence from B (0,ml,m^ ) is canceled by the divergence from .B (0, m 2 , m~ ), 
making finite the contribution from bottom-sbottom loops to this sub-matrix, as it 
should be, since the mass is fully "calculable" . 

The second most important contribution to the radiatively induced neutrino 
mass usually comes from charged-scalar/charged-fermion loops 26]. Since all possible 
topologies of Fig. (3) contribute to this loop the structure of the contribution from 
charged Higgs/slepton loops is more complex than that of the bottom-sbottom loop 
considered above. However, the same topology as for the sbottom/bottom loop also 
contributes to the charged scalar loop. It leads to a final expression similar to eq. 
(|5ip . with appropriate replacements, which is good enough for an order-of-magnitude 
estimate of the charged scalar loop. 

4-1. Results for the solar mass scale 

We give a discussion of the analytical versus numerical results of the solar mass scale 
first. In Fig. (jSJ) we show a comparison of approximate and exact calculation for two 
different numerical data sets. In both figures we show the ratio of the approximate-over- 
exact solar neutrino mass parameter m^ pr / m e ^ act versus Am s 2 OL in eV 2 , where m^f pr is 
the approximate loop calculation involving the bottom-sbottom and the charged scalar 
loop, while m e J^ act is the exact numerical computation taking into account all loops. The 
set to the left called "Ntrl" contains neutralinos being the LSP, while the set to the right 
(Stau) has the charged scalar tau as LSP. 
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Am 2 „ n . [eV 2 ] Am 2 ,™ [ ey2 ] 



Figure 5. (m^ pr /m^ act ) versus Amf 0L [eV 2 ] for the set Ntrl (left) and the set 
Stau (right). ra^ pr is the sum of the bottom-sbottom and charged scalar loops, while 
m exact j s |_ ne numer j ca i result for all loops. In case of LMA the approximation works 
always better than 10 %. 



We have found numerically that the terms proportional to e» x ij in the self energies 
in Eq. (|51|) give the most important contribution to m V2 in the bottom-sbottom loop 
calculation in most points of our sets. If these terms are dominant one can find a very 
simple approximation for the bottom-sbottom loop contribution to m U2 . It is given by 

m U2 ~ -jL sin(29~ b )m b ABl^ M + ® . (55) 

Eq. ()55)l works surprisingly well for almost all points in our data sets. 

The more complicated structure of the charged scalar loop makes it difficult to give 
a simple equation for m V2 similar to Eq. ()55)1 for the bottom-sbottom loop. However, 
we note that Eq. 1)5 5 j). with appropriate replacements, allows us to estimate the typical 
contributions to the charged scalar loop within a factor of ~ 3. However, such an 
estimate will be biased toward too small or too large m„ 2 depending mainly on which 
SUSY particle is the LSP Q- 



4-2. Analytical approximation for the solar mixing angle 

In the basis where the tree-level neutrino mass matrix is diagonal the mass matrix at 
one-loop level can be written as 

\ 



m v = K ( ° )T ^K (0) 



cieiei cieie 2 Cieie 3 
'12 



+ • ■ ■ (56) 



Ci6 3 6i cie 3 e 2 c |A| 2 + cie 3 e 3 J 

where the ?j were defined before in Eq. (jHHJl . Coefficients Co and c\ contain couplings 
and supersymmetric masses. Since they cancel in the final expression for the angle their 
exact definition is not necessary in the following. Dots stand for other terms which 
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we will assume to be less important in the following. This matrix can be diagonalized 
approximately taking in account that 



x 



c |A| 2 



< 1 



(57) 



Then 



x 



c |Ap 



^1 



if 



X ■ 



x ■ 



x ■ 



ei£2 

if 
£262 



:r ■ 



X ■ 



if 

6362 



;r ■ 



£i£3 \ 

if 
£263 



1 + 

|e|* 



(58) 



The rotation matrix that diagonalizes m„ in Eq. (|58|) can be written as 
Vjm v V v = diag(mi, m 2 , m 3 ) 



where 



V 1 



ei,i 
e 2 ,i 
\ e 3,i 

The lepton mixing matrix is then given by 

T 



ei )2 
e2,2 
e 3 ,2 



e l,3 

e2,3 

e 3,3 / 



1/ 



The expression for the solar mixing angle can be obtained from: 
tan 2 6>c 



'SOL 



u 2 
u 2 



(59) 



(60) 



(61) 



(62) 



From the above equations we obtain the very simple expression for the solar mixing 
angle, 



tan 2 6U, = -Jr 



r2 



(63) 



This formula is a very good approximation if the one-loop matrix has the structure 
€i x €j, as is the case of the bottom-sbottom loop if m U3 3> m U2 , as illustrated in Fig. El 
In the left panel we show a calculation comparing for all points the approximate 
to the exact solar angle in the set with neutralino LSP, while the right panel shows a 
subset of points with the cut sin(26'j ) )A£>Q 2ri > 0.02. Note that this cut is designed so 
as to favor points in which there is a sizeable bottom-sbottom loop contribution to the 
full one-loop neutrino mass. One sees from the right panel that for this case the true 
solar angle is well approximated by our analytical formula. Note finally that eq. (jo^|) 



will fail completely, if A M = A T and = 
origin of the "sign condition" discussed in |2f 



since then e| 



0, see Eq. (JHSJ)- This is the 
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« 5 

CO 



< o 




0.2 0.4 0.6 0.8 1 

2 -Exact 

tan SOL 



Figure 6. (tan 2 



r , o A exacts . 9 /i 

/tan v SOL ) versus tan st 



On the left panel the 



darker region contains over 90% of the points in our sample. In the right panel the 
points in the region shown satisfy the cut sin(26'j ) )Ai3Q 2ri > 0.02 . 



5. Testing neutrino properties at high energy accelerators 

Since R-parity is broken in our model, the lightest supersymmetric particle is unstable 
and decays. This leads to the exciting possiblity to test the bilinear model at future 
colliders, such as the LHC or a possible Linear Collider. 



The principle idea of such a collider test 39, 41, 43] is easily understood: Bilinear Im- 



parity breaking leads to mixing between particles and sparticles with the same quantum 
numbers, as discussed above extensively for the case of neutrinos/neutralinos. This 
mixing, however, is not arbitrarily different for each particle/sparticle species. In fact, 
the bilinear model has just six new parameters, which we choose to be e« and A i; 
compared to the MSSM. Essentially five of these six can be fixed from neutrino physics. 

Thus, if the MSSM parameters were known, all mixing effects could be calculated 
and thus all decay properties of the LSP would be fixed - apart from the effects of the 
last unknown parameter. In reality, however, the MSSM soft SUSY breaking parameters 
are completely unknown. The approach taken in 39, 41, therefore is to calculate 



ratios of branching ratios of different decays. By taking ratios one essentially scales 
out the unknown MSSM parameters approximately and obtains observables which are 
proportional to either Aj/Aj or 6j/e 3 - (or some weird combination thereof). Which ratio 
one measures depends of course on the final state and the LSP under consideration. 
Since ratios of Aj's (or e^'s) are correlated with the neutrino angles, as discussed above, 
fixing neutrino angles from experimental data therefore gives definite predictions for 
some ratios of branching ratios. 

One example is shown in Fig. where ratio of branching ratios of neutralino LSP 
decays are plotted. Note that Br(fiqq') /Bi^rqq') is directly proportional to ta.n 2 (9 atm ), 
i.e. should be near ~ 1 according to current neutrino data. The spread in the points 
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Figure 7. Ratio of branching ratios of neutralino LSP decays. To the left: 
Bi(/j,qq') /Br(rqq') as a function of tan 2 ((9 atm ). To the right: Br(eqq') /~Br(rqq') as 
a function of U^ s . Whether bilinear R-parity breaking SUSY is responsible for 
atmospheric neutrino oscillations can be checked by such a measurement. Note that 
the spread of the points is entirely due to the unknown MSSM parameters. Even a 
moderately accurate input of MSSM parameters will lead to much sharper predictions 
for these decays. 



is due to the unknown MSSM parameters. Of course, once SUSY is discovered these 
unknowns will be measured allowing for much sharper tests of the model than indicated 
in Fig. [7| 

A second example is shown in Fig. where we show Br(fi — > *Y^ ve )l Br(fi — > 
*Y^v\i) as a function of (ei/e2) 2 (left) and as a function of tan 2 (# ) (right). Obviously 
this ratio is strongly correlated with the solar angle and thus, if scalar taus turn out to 
be the LSP, such a measurement would provide an excellent check of the bilinear model 
as the origin of the solar neutrino mass scale. 

With the LSP unstable, in principle any spart.de can be the LSP. In the 
remaining candidates have been discussed: Charginos, scalar quarks, gluinos and scalar 
neutrinos. The main conclusion of ji^ is that whichever SUSY particle is the LSP, 
measurements of branching ratios at future accelerators will provide a definite test of 
bilinear R-parity breaking as the model of neutrino mass. We just mention that chargino 
LSPs would be more sensitive to atmospheric neutrino physics (as are neutralinos) while 
the other LSP candidates mentioned above show more dependence on the solar neutrino 
angle. 

6. Discussion and conclusions 

We have presented a brief review of the idea that supersymmetry with explicit bilinear 
breaking of R-parity is the origin of neutrino masses and lepton mixing. The bilinear R- 
parity breaking (BRpV) model is the simplest extension of the minimal supersymmetric 
standard model (MSSM) which includes lepton number violation. We have seen how 
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Figure 8. Ratio of branching ratios of scalar tau LSP decays. To the left: 
Br(fi — > ve ) /Br(fi — > yj v iA as a function of {tijtri) 1 ■ To the right: as a function 
of tan 2 (00 ). Whether bilinear R-parity breaking SUSY is responsible for solar neutrino 
oscillations can be checked by such a measurement if scalar taus arc the LSP. 

it leads to a successful phenomenological model for neutrino oscillations, in accordance 
to present neutrino data. The pattern of neutrino masses is hierarchical, with the 
atmospheric mass scale arising at the tree level whereas the solar scale is induced from 
calculable loop corrections. We saw how, in contrast to seesaw models, the BRpV model 
can be probed at future collider experiments, like the LHC or the NLC Indeed we have 
discussed how, irrespective of the supersymmetric particle which is the lightest, its decay 
pattern will be directly related with the lepton mixing angles determined in low energy 
neutrino experiments. 
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